Theoretical interpretation of the experimental electronic structure of lens shaped, 

self-assembled InAs/GaAs quantum dots 
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We adopt an atomistic pseudopotential description of the electronic structure of self-assembled, 
lens shaped InAs quantum dots within the "linear combination of bulk bands" method. We 
present a detailed comparison with experiment, including quantities such as the single particle 
electron and hole energy level spacings, the excitonic band gap, the electron-electron, hole-hole 
and electron hole Coulomb energies and the optical polarization anisotropy. We find a generally 
good agreement, which is improved even further for a dot composition where some Ga has 
diffused into the dots. 
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I. INTRODUCTION: USING THEORY AS A 
BRIDGE BETWEEN THE STRUCTURE AND 
THE ELECTRONIC PROPERTIES OF 
QUANTUM DOTS 

Self-assembled, Stranski-Krastanow (SK) grown semi- 
conductor quantum dots have recently received consider- 
able attention as they exhibit a rich spectrum of phenom- 
ena including quantum-confinement^], |^, 




exchange- 




Coulomb charging/blockadej 
and multi-exciton transitions! 
the past few years a considerable number of high quality 
measurements of the electronic level stucture of these dot 
systems have been performed, using photoluminescence 
(PL) go] If, H, U, |l|, p hotoluminescence lu- 

minescence excitation (PLE)|4, |14|, capacitance) 
HI and far infra red (FIR) spectroscopy]^, |l], p2|, 



24|. 

These measurements have been able to determine the 
electronic level structure to relatively high precision. In 
parallel with these measurements, several groups have 
also attempted to measure th e g eometry and composi- 
tion of these dotsjl^, [l6[ 




Eq, 27 . So far, however, 



these measurements have failed to provide details of the 
shape, size, inhomogeneous strain and alloying profiles 
to a similar level of accuracy to that in which the elec- 
tronic structure has been determined. As a result, the 
size of the dots were often used as adjustable parame- 
ters in models that fit experimental spectra. For exam- 
ple, using a single-band effective mass model, Dekel et. 
q/p^] defined an "effective shape" (cuboid) and "effec- 
tive dimension" that reproduced the measured excitonic 
transitions. Similar "parabolic dot" models have been 
assumed by Hawrylak et. a/[0. 

The accuracy of single-band and multi-band effective 
mass methods was recently examined in a series of pa- 
,§| §§. 
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In these works, the shape, size 
and composition of nanostructures were arbitrarily fixed, 
and the electronic structure was evaluated by succes- 
sively improving the basis set, starting from single-band 
methods (effective-mass), going to six and eight band 
methods (k.p), and finally, using a converged, multi- 



band approach (plane- wave pseudopotentials) . It was 
found that conventional effective-mass and k.p methods 
can sometimes significantly misrepresent the fully con- 
verged results even when the shape, size and composi- 
tion was given. The observed discrepancies were both 
quantitative (such as band gap values, level spacings, 
Coulomb energies) and qualitative (absence of polariza- 
tion anisotropy in square based pyramidal dots [[32], miss- 
ing energy levels pp|| ) . As a result of these limitations 
these methods may not offer a reliable bridge between 
the electronic structure and atomic structure. 

In this paper, we offer a bridge between recent mea- 
surements of the electronic structure and measurements 
of the atomic structure of the dots using accurate the- 
oretical modeling. Modeling can determine if the cal- 
culated electronic structure resulting from an assumed 
shape, size, strain and alloying profiles agrees with the 
measured electronic structure or not. A theory that can 
perform such a "bridging function" must be accurate and 
reliable. The pseudopotential approach to this problem 
qualifies, in that any discrepancy between the predicted 
and measured electronic properties can be attributed to 
incorrectly assumed shape, size or alloying profile. We 
have studied a range of shapes, sizes and alloy profiles 
and find that a lens-shaped InAs dot with an inhomo- 
geneous Ga alloying profile is in closest agreement with 
current measurements. In the following sections we at- 
tempt to provide a consistent theoretical interpretation 
of numerous spectroscopic properties of InAs/GaAs dots. 



II. OUTLINE OF THE METHOD OF 
CALCULATION 

We aim to calculate the energy associated with vari- 
ous electronic excitations in InAs/GaAs quantum dots. 
These energies can be expressed as total energy differ- 
ences and require four stages of calculation: 

(i) Assume the shape, size and composition and com- 
pute the equilibrium displacements: We first construct 
a supercell containing both the quantum dot and sur- 
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rounding GaAs barrier material. The shape, size and 
composition profile are taken as input and subsequently 
refined. Sufficient GaAs barrier is used, so that when 
periodic boundary conditions are applied to the sys- 
tem, the electronic and strain interactions between dots 
in neighboring cells is negligible. The atomic posi- 
tions within the supercell are then relaxed by minimiz- 
ing the strain energy described by an atomistic force 
field |33|, 34| including bond bending, bond stretching 
and bond bending-bond stretching interactions (see sec- 
tion [iIIA). An atomic force field is similar to continuum 
elasticity approaches |54j in that both methods are based 
on the elastic constants, {CV,}, of the underlying bulk 
materials. However, atomistic approaches are superior to 
continuum methods in two ways, (a) they can contain an- 
harmonic effects, and (b) they capture the correct point 
group symmetry, e.g. the point group symmetry of a 
square based, zinc blende pyramidal dot is Civ-, since the 
[110] and [110] directions are inequivalent while contin- 
uum methods ]34|, find C^ v . M ore de tails of the atomistic 
relaxation are given in section III A . 

(ii) Setup and solve the pseudopotential single-particle 
equation: A single-particle Schrodinger equation is set up 
at the relaxed atomic positions, {R nQ } 



<-§v> 



(r - R„ Q )}V> 4 (r) = a ipi(r) 



(1) 



The potential for the system is written as a sum of 
strain-dependent, screened atomic pseudopotcntials, v a , 
that are fit to bulk properties extracted from expe riment 
and first-principles calculations (see section III B|) . The 



Schrodinger equation is solved by expanding ip in a lin- 
ear combination of bulk states, (f> n k, from bands, n, and 
k-points, k, 



</>i(r,e) 



n,k 



n.k 



<?!>„fc(r,e) 



(2) 



taken at a few strain values. The solution of Eqs. (|lj) 
and ^ provides the level structure and dipole transi- 
tion matrix elements. More details on the solution of the 



Schrodinger equation are given in section IIIC. 

(iii) Calculate the screened, inter-particle many-body 
interactions: The calculated single particle wavefunc- 
tions are used to compute the electron-electron, electron- 
hole and hole-hole direct , J ee , J e h, Jhh, an d exc hange 
K ee , K e h, Khh Coulomb energies (see section HID). 

(iv) Calculate excitation energies as differences in to- 
tal, many particle energies: For example, the difference 
between the total energy En [/ijej] of a dot with a hole in 
level h and an electron in level eo and the total energy 
-Eoo[^o e ol of the unexcited dot is 



Eoo[h^e^} 



(e eo - e ho ) 

Je h + 2K eoho 5sfi 



(3) 



where (in the absence of spin-orbit coupling) 63,0 = 1 
for triplet states, and for singlet states. Analagous 



express ions e xist for electron-addition experiments (see 
section III D ). 

The main approximations involved in our method are: 

(a) the fit of the pseudopotential to the experimental 
data of bulk materials is never perfect (see Table ||) and 

(b) we neglect self-consistent iterations in that we assume 
that the screened pseudopotential drawn from a bulk cal- 
culation is appropriate for the dot. Our numerical con- 
vergence parameters are (i) the size of the GaAs barrier 
separating periodic images of the dots, and (ii) the num- 
ber of bulk wavefunctions used in the LCBB expansion 
of the wavefunctions. To examine the effects of these 
approximations and convergences on the ultimate level 
of accuracy that can be obtained with our methodology 
we have first applied the se me thods to an InGaAs/GaAs 
quantum well(see section III E ) , where experimental mea- 
surements of the shape, size, compositio n and transition 
energies are more established (see section III E). We next 
describe the details of our method. 



TABLE I: Fitted bulk electronic properties for GaAs 
and InAs using the screened atomic pseudopotentials, 
in Eq.(^). The hydrostatic deformation potential of the 
band gap and Ti§ v levels are denoted by a gap and ar 15t , . 
The biaxial deformation potential is denoted by b and 
the spin-orbit splittings at the Ti^ v and L\ v points are 
denoted by Aq and Ai. 



Property 


EPM 


GaAs 

ExptjiJ 


EPM 


InAs 

Exptfiil 


Egap 


1.527 


1.52 


0.424 


0.42 




-2.697 


-2.96 


-2.330 


-2.40 


Ex la 


1.981 


1.98 


2.205 


2.34 


Ex 3a 


2.52 


2.50 


2.719 


2.54 




-1.01 


-1.30 


-5.76 


-6.30 


E Llc 


2.36 


1.81 


1.668 


1.71 


* 


0.066 


0.067 


0.024 


0.023 


m* hh [W0] 


0.342 


0.40 


0.385 


0.35 


m* hh [lll] 


0.866 


0.57 


0.994 


0.85 


m,\[100] 


0.093 


0.082 


0.030 


0.026 


@>gap 


-7.88 


-8.33 


-6.79 


-5.7 


"r 15 „ 


-1.11 


-1.0 


-0.826 


-1.0 


b 


-1.559 


-1.7 


-1.62 


-1.7 


Ac 


0.34 


0.34 


0.36 


0.38 


Ai 


0.177 


0.22 


0.26 


0.27 



III. DETAILS OF THE METHOD OF 
CALCULATIONS 

A Calculation of equilibrium atomic positions for a 
given shape 

To calculate the relaxed atomic positions within the 
supercell, we use a generalization (G-VFF) of the origi- 
nal valence force field (VFF)|33| model. Our implemen- 
tation of the VFF includes bond stretching, bond angle 
bending and bond- length/bond- angle interaction terms 
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in the VFF Hamiltonian. This enables us to accurately 
reproduce the Cn, Cn and C 44 elastic constants in a 
zincblende bulk material. We have also included higher 
order bond stretching terms, which lead to the correct 
dependence of the Young's modulus with pressure. The 
G-VFF total energy can be expressed as: 



.(2) 



Ad? 



££ 



3/? 



k>j 



9° u , 

jik ij 



nn 



i k>j 

cos6% k dld° k } , 



where Ag? 2 ,- 



[(Ri - Rj 



d ' 



]/d\ 



(4) 



Here R, is the 



coordinate of atom i and df-j is the ideal (unrelaxed) bond 
distance between atom types of i and j. Also, 8j ik is 
the ideal (unrelaxed) angle of the bond angle j — i — k. 
The YJ n ™ 1 denotes summation over the nearest neigh- 
bors of atom i. The bond stretching, bond angle bend- 
ing, and bond- length/bond-angle interaction coefficients 

otij{= a), 'jik, a jik are related to the elastic constants 
in a pure zincblende structure in the following way, 

Cu +2C 12 = i-{3a + [3-6a) 

C\x — Cl2 = — —p 
do 

_ \/3 [{a + P)(a0 -a 2 )-2a 3 + 2al3a} 
° 44 _ ~do~ (a + /3 + 2a) 2 (&) 

The second-order bond stretching coefficient is re- 
lated to the pressure derivative of the Young's modulus 
by where B — (Cn + 2Ci 2 )/3 is the Young's mod- 
ulus. Note that in the standard [Q VFF which we have 
used previously^, [56], [jTjthe last terms of Eq.([|) are 
missing, so a = in Eq.(fcl). Thus there were only two 
free parameters (a, (3) and therefore three elastic con- 
stants could not, in general, be fit exactly. The G-VFF 
parameters and the resulting elastic constants are shown 
in Table || for GaAs and InAs crystals. For an InGaAs al- 
loy system, the bond angle and bond- length/bond-angle 
interaction parameters f3, a for the mixed cation Ga-As- 
In bond-angle are taken as the algebraic average of the 
In-As-In and Ga-As-Ga values. The ideal bond angle 
6j ik is 109° for the pure zincblende crystal. However, to 
satisfy Vegas's law for the alloy volume, we find that it 
is necessary to use 0Q a _ As _ In — 110.5° for the cation 
mixed bond angle. 

As a simple test of this G-VFF for alloy systems, 
we compared the relaxed atomic positions from G- 
VFF with pseudopotential LDA results for a (100) 



TABLE II: Input G-VFF parameters a, (3, a to Eq.(@) 
and their resulting elastic constants Cn, Cn and C44 
(10 3 dyne/cm) 





a 





a 


a (2) 


Cn 


C12 


C44 


GaAs 


32.153 


9.370 


-4.099 


-105. 


12.11 


5.48 


6.04 


InAs 


21.674 


5.760 


-5.753 


-112. 


8.33 


4.53 


3.80 



(GaAs) 1/ (InAs) 1 superlattice where the c/a ratio is fixed 
to 1, but we allow energy minimizing changes in the over- 
all lattice constant (a eq ) and the atomic internal degrees 
of freedom (u eq ). We find a^ q DA = 5.8612 A and u^ q DA = 

5.8611 

0.2305. In comparison the original 



0.2305, while the G-VFF results are af q - VFF 
A and u?~ VFF = 



VFF yields a v eq ' 



VFF 



5.8476 A and u v e FF 



0.2303. 



B The Empirical Pseudopotential Hamiltonian 



We set up the single-particle Hamiltonian as 

i? = -fv 2 +£v Q (r-R„ a ) , (6) 



where R„ Q is the G-VFF relaxed position of the n th atom 
of type a. Here v a (r) is a screened empirical pseudopo- 
tential for atomic type a. It contains a local part and a 
nonlocal, spin-orbit interaction part. 

The local potential part is designed to include depen- 
dence on the local hydrostatic strain Tr(e): 



^ c (r; £ )=tC(r;0)[l+ 7Q Tr( e )] 



(7) 



where the 7 a is a fitting parameter. The zero strain po- 
tential V^(r; 0) is expressed in reciprocal space q as 



= a o(? 2 - a l)/[ a 2< 



(8) 



The local hydrostatic strain Tr(e) for a given atom at 
R is defined as £Ir/Qq — 1, where £Ir is the volume of 
the tetrahedron formed by the four atoms bonded to the 
atom at R. £Iq is the volume of that tetrahedron in the 
unstrained condition. The need for explicit dependence 
of the atomic pseudopotential on strain in Eq. (ffl) results 
from the following: While the description in Eq. (p) of the 
total pseudopotential as a superposition of atomic poten- 
tials situated at specific sites, {R na }, does capture the 
correct local symmetries in the system, the absence of a 
self-consistent treatment of the Schrodinger equation de- 
prives the potential from changing in response to strain. 
In the absence of a strain-dependent term, the volume 
dependence of the energy of the bulk valence band maxi- 
mum is incorrect. While self-consistent descriptions show 
that the volume deformation potential a v = dE v /dhiQ of 
the valence band maximum is negative for GaAs, GaSb, 
InAs, InSb and for all II- VI this qualitative behavior can 
not be obtained by a non-self-consistent calculation that 
lacks a strain dependent pseudopotential. 
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The nonlocal part of the potential describes the spin- 
orbit interaction, 

H so = J2V«°(Rn a ) (9) 

not 

= EE v i s ,> - . s(i\ Rna , 

na I 

where \l)n na is a projector of angular momentum I cen- 
tered at R na , L is the spatial angular momentum opera- 
tor, S is the Dirac spin operator, and Vf^{r) is a potential 
describing the spin-orbit interaction. 

In Eq(||), the kinetic energy of the electrons has been 
scaled by a factor of (3. The origin of this term is as 
follows: In an accurate description of the crystal band 
structure, such as the GW methodQ, a general, spa- 
tially non-local potential, V(r,r'), is needed to describe 
the self-energy term. In the absence of such a term the 
occupied band width of an inhomogcncous electron gas 
is too large compared to the exact many-body result. 
To a first approximation, however, the leading effects of 
this non-local potential, V(r, r'), can be represented by 
scaling the kinetic energy. This can be seen by Fourier 
transforming V(r, r') in reciprocal space, q, then making 
a Taylor expansion of q about zero. We find that the 
introduction of such a kinetic energy scaling, (3 permits a 
simultaneous fit of both the effective masses and energy 
gaps. In this study, we fit j3 — 1.23 for both GaAs and 
In As. 

The pseudopotential parameters in Eqs(Q) and (||) were 
fitted to the bulk band structures, experimental deforma- 
tion potentials and effective masses and first-principles 
calculations of the valence band offsets of of GaAs and 
InAs. The alloy bowing parameter for the GalnAs band 
gap (0.6 eV) is also fitted. The pseudopotential param- 
eters are given in Table 111 and their fitted properties 
are given in Table ||3£j. We sec that unlike the LDA, 
here we accurately reproduce the bulk band gaps and 
the bulk effective masses. One significant difference in 
our parameter set, to that used in conventional k.p stud- 
ies, is our choice of a negative magnitude for the valence 
band deformation potential, a v , which we have obtained 
from LAPW calculations Bol . 



TABLE III: Parameters for the GaAs and InAs screened 
atomic pseudopotentials, in Eq.(Q). This potential re- 
quires a plane wave cutoff of 5 Ryd. 



Parameter 


In 


Ga 


As (InAs) 


As (GaAs) 


ao 


644.13 


432960 


26.468 


10.933 


ai 


1.5126 


1.7842 


3.0313 


3.0905 


a.2 


15.201 


18880 


1.2464 


1.1040 


a3 


0.35374 


0.20810 


0.42129 


0.23304 


&4 


2.1821 


2.5639 


0.0 


0.0 



The present InAs and GaAs pseudopotentials have 
been systematically improved relative to our previous 
InAs and GaAs potentials ||, |f| |l|, ||, ||, although 
the functional form has remained the same. Firstly, the 



pseudopotentials for InAs and GaAs used in Ref.Kq, 42] 
did not include the spin-orbit interaction. In Rcfs.^Tj. 3(j. 
|]| we used potentials that included the spin-orbit inter- 
action, but were not able to simultaneously, accurately 
fit the electron effective and the zone center band gap, 
due to the lack of the above (3 parameter. The potential 
used here is identical to that used in Refs.M. B2|. 



C Calculating the single particle eigenstates 

One could use a straight forward expansion of the sin- 
gle particle wavefunctions in a plane wave basis set, as 
we have previously done in Refs.[|^, |3(| [35|. However, as 
was shown in Refs.j|| H @ ; 

a more economical repre- 
sentation is to use the Linear Combination of Bulk Bands 
(LCBB) method||, ||, @. Within the LCBB the eigen- 
states of the pseudopotential Hamiltonian are expanded 
in a basis of bulk Bloch orbitals 



V'i(r) = EE c ^,fc Us ' n > k 



(10) 



n ,k 



where u SiIlj fc(r) is the cell periodic part of the bulk Bloch 
wavefunction for structure, s, at the n th band and the 
k th k-point, k. These states form a physically more in- 
tuitive basis than traditional plane waves therefore the 
number of bands and k-points can be significantly re- 
duced to keep only the physically important bands and 
k-points (around the T point in this case). This method 
was recently generalized to strained semiconductor het- 
erostructure systems^] and to include to spin-orbit 
interaction j43). In this paper use an LCBB basis de- 
rived from four structures, s. These structures are (i) 
unstrained, bulk InAs at zero pressure, (ii) unstrained, 
bulk GaAs at zero pressure, (iii) bulk InAs subjected to 
the strain value in the center of the InAs dot, and (iv) 
bulk InAs subjected to the strain value at the tip of the 
InAs dot. By interpolating the strain profile between 
these four structures, the basis is able to accurately de- 
scribe all the strain in the system. The wavevectors, {k}, 
used here include all allowed values within 167r/L of the 
zone center, where L is the supercell size. For calcula- 
tions of electron states, the band indices, n, include only 
the band around the Ti c point. For the hole states we 
also include the three bands around the Ti^ v point. This 
basis set produces single particle energies that are con- 
verged with respect to basis size, to within 1 meV. 



D Constructing the energies of different electronic 
configurations 

Using screened Hartree Fock theory, the energy asso- 
ciated with loading N electrons into a quantum dot can 
be expressed p5| as 



En = + s 



pol 



(11) 



5 



where q are the single-particle energies and E? oZ are the 
polarization self-energies of the i th electron state , J^ e 
and Kfj are the direct and exchange Coulomb integrals 
between the i th and j th electronic states and rii are the 
occupation numbers (J2i n-i — N). As shown in Ref.Q, 
for free standing, colloidal quantum dots the dielectric 
constant inside the dot is dramatically different to that 
outside (vacuum) and hence the polarization self-energy, 
, is very significant (~1 eV). For self assembled InAs 
dots embedded in GaAs, the dielectric constants of InAs 
and GaAs are similar (£ m = 12.3, 10.6) and we calculate 
this term as ~1 meV, and hence we choose to neglect 
it here. The direct and exchange Coulomb energies, are 
defined as 



J i 



\Uri)\ 2 l^(r 2 )| 2 
e(ri - r 2 )|ri - r 2 



dr 1 dr2 



V>*(n) y>,(r 2 ) y>*(r 2 ) y>j(n 

e(ri - r 2 )|ri - r 2 | 



-dvidv 2 (12) 



where e is a phenomenological, screened dielectric 
function ]4l|] containing a Thomas Fermi electronic com- 
ponent and an ionic component from Ref . Jill . Our ex- 
change automatically includes both short and long range 
components. 

Denoting electron levels as eo,ei,e 2 ..., hole levels as 
ho, hi, hi-., and the number of electrons and holes as N 
and M, the total energy, Emn, is 



E M N = J2 ~ e hi m i + YP^O 1 - K ^) m * m 3 
i i<j 

i i<j 

, (13) 



where n% and rrii are the electron and hole occupation 
numbers respectively and J2i Ui = N and ^ m, = M. 
Using Eq.(fl3"|), in the strong confinement regime where 
kinetic energy effects dominate over the effects of ex- 
change and correlation, an exciton involving electrons ex- 
cited from hole state i to electron state j can be expressed 
as 



jpexciton 



(14) 



To study charged dots, if one assumes the electron states 
are occupied in order of increasing energy (Aufbau prin- 
ciple), the total energy of a dot charged with TV electrons, 
Eqn, is 



oo L e o 

bi H] 

02 [eo 
^03 [e 2 e\] 
Eqa [ele\] 



]=0 

] = £ e 

] = 2e eo + Je ,e 

(2^eo "T" Ee% ) "T" [^eo,eo ~t~ 2</ e0jei ] K e 



(2e eo 
2K en 



2e„ 



[Je, 



Je 



4 J, 



(15) 



As indicated in section |l| our wavefunctions, {ipi}, are 
not iterated to self-consistency. This affects the magni- 
tude of the direct and exchange Coulomb integrals. We 
have previously examined the accuracy of this perturba- 
tive treatment for colloidal InAs dots by comparing the 
non-self-consistent Coulomb energy with that obtained 
self consistently [47 . The differences were negligible. 



E Quantum well tests 

To test the above methods, we first calculated the 
energy levels in a quantum well, and compared the re- 
sults with experiment. In Fig. |l|(a), we compare the 
calculated electron-heavy hole transition energies for a 
96 A Ino. 2 4Gao.76As quantum well inside a GaAs matrix. 
The peaks in the experimental spectra occur|Q at 1.275, 
1.395 and 1.538 eV. Our calculated transitions occur at 
1.290, 1.404 and 1.545 eV respectively. Figure 0(b) com- 
pares the band gap of a Ino. 22 Gao.7sAs quantum well as 
a function of its thickness. The measured band gaps|49| 
for quantum wells with thicknesses of 6 and 18 ML are 
1.458 and 1.351 eV. Our calculated values are 1.466 and 
1.366 eV. 



IV. PHYSICAL QUANTITIES TO COMPARE 
WITH EXPERIMENT 



The quantities we use to characterize the electronic 
structure are illustrated in Fig. |^ which shows a 
schematic layout of the electron and hole single-particle 
energy levels in a quantum dot. Assuming that all levels 
are spatially nondegenerate (thus having only spin degen- 
eracy), we mark the electron levels as eo, e\, e 2 .... and the 
hole levels as ho,hi,h 2 . The level eo is sometimes called 
"s-like", whereas ei and e 2 are called "p-like", and e% 
and e4 are called "d-like" . Since the GaAs environment 
of the InAs dots is largely unstrained, it is convenient to 
set as a reference energy the VBM of GaAs as E = 0, 
and the CBM of GaAs as E=1520 meV. All energy levels 
can be referenced with respect to these band edges. 
For the electron levels, the quantities that we consider 
are: 

(i) The number of dot-confined electron states, N e . 

(ii) The spacing S sp = e ei — e eo between "s-like" and 
"p-like" electron states. 

(iii) The splitting S pp = e e2 — e ei between the "p-like" 
electron states. 

(iv) The spacing S pc i — e E3 — e e2 between "p-like" and 
"d-like" electron states. 

(v) The "binding energy" of the first electron level, 
eo, with respect to the GaAs conduction band minimum, 

AE(e) = EcaAs,CBM — Ceo- 

(vi) The position of the bottom of the band for the 
2D InAs "wetting layer" (WL) with respect to the GaAs 



CBM, AE$ L = E GaAStC BM 



E 



(e) 
WL- 



Fig. 1 Williamson ef. al. 
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(a) QW Excitonic Transitions 
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FIG. 1: (a) A comparison of EPM calculated and mea- 
sured electron to heavy hole transition energies in a 96 
A Ino.24Gao.76As quantum well embedded inside a GaAs 
matrix. The vertical lines mark the positions of the EPM 
calculated transitions, (b) The calculated band gap of an 
Ino.22Gao.78As quantum well as a function of its thick- 
ness. 



(vii) Inter-electron direct J^ ej and exchange K^f ej 
Coulomb energies. 
For the hole levels we consider are: 

(i) The number, Nh, of dot-confined hole states. 

(ii) The intra-band spacings of the hole levels, 5^ = 

thj — Chi ■ 

(iii) The "binding energy" of the first hole level, 
ho, with respect to the GaAs valence band maximum, 
—AE(h) — Ego,As,vbm + e h - 

(iv) The position of the top of the band for the 2D InAs 
"wetting layer" (WL) with respect to the GaAs VBM, 

^ e wl = —Ego,As,vbm + e wl- 

(v) Inter-hole direct J^ hj and exchange K^ h , 

Coulomb energies. 

Finally, for the recombination of electrons and holes, we 
consider: 



(i) The excitonic energies, Eff clton , as defined in 
Eq.(|14|). By subtracting calculated values for the single 
particle energies e Ej and from measured optical exci- 
tation energies one can estimate the electron-hole direct 
Coulomb energies Jh^j- 

(ii) The ratio of absorption intensities for light polar- 
ized along [110] and [110] directions, defined as 



A = 



-P[110] _ < ^e \r[l W \\4>h 

0^ 



-P[110] < V'eol^lTO]!^ 



(16) 



This ratio can deviate from unity due to three reasons; 
(a) The dots has different dimensions in the [110] and 
[110] directions. We refer to this as the the "geometric 
factor" . (b) The atomistic zincblende symmetry makes 
the two directions symmetry inequivalent even if the 
lengths along the two directions are equal. We refer 
to this as the "atomic symmetry factor". One mani- 
festation of this affect is that if the strain is calculated 
atomistically, it is different in the two directions even in 
the absence of a geometric factor Q. (c) A piezoelec- 
tric field that breaks the symmetry. Previous studies [T^] 
have shown that this effect is negligible in InAs/GaAs 
dots so we will neglect it here, k.p calculations neglect 
the "atomic symmetry" factor (except for the small effect 
of strain asymmetry), but retain the "geometric factor". 
Pseudopotential calculations retain both effects. For ex- 
ample, in a square based pyramid (where by definition 
the "geometric factor" does not contribute), k.p produces 
A = 1 while pseudopotential theory gives A = 1.2 (see 
Table |V|) . This shows that there is not a simple mapping 
from dot shape to polarization anisotropy, A. 

(iii) Excitonic dipole: As the center of the electron 
and hole wavefunctions do not exactly coincide with each 
other, it is possible that an exciton will exhibit a de- 
tectable dipole moment, 



dh z ,e 3 = (iphil^hi) ~ {^ ej \f\tp ej ) 



(17) 



The quantities defined above characterize the elec- 
tronic structure. Next, in section [v|, we will next provide 
all of these quantities from our calculations, and then 
in section [V| we will extract measured values of these 
quantities from the available experiments. 



V. THEORETICAL RESULTS 

The electronic structure of a series of GalnAs/GaAs 
self- assembled quantum dots was calculated using the 
methodology described in Section |l[ We have chosen 
to focus on the well established "lens shaped" dot geom- 
etry from Refs.§, |, @, @H |H @. The shape of 
this dot is shown in Fig. pT The profile is obtained by 
selecting the section of a pure InAs sphere that yields a 
circular base with diameter 252A and a height of 35 A. 
The main experimental uncertainty about this dot is the 
composition profile. It is not known if the dots are pure 



Fig. 2 Williamson ef. a/. 



Fig. 3 Williamson ef. a/. 



GaAs 



GaJn^As 



GaAs 



AWL(e) 




FIG. 2: Schematic single particle electron and hole en- 
ergy levels for lens shaped, InAs quantum dots embedded 
within GaAs. 
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FIG. 3: Assumed model geometry of the lens and pyra- 
midal shaped quantum dots. 



InAs or if Ga has diffused into the dots. For compari- 
son, we also calculate the electronic structure of a square 
based InAs pyramid with a base of 113A and a height of 
56A. This is not believed to be a realistic geometry, how- 
ever, it has been used as a benchmark for many previous 
theoretical calculations |32| |||, |36|, [30[ |5|] and we include 
it here for comparison purposes. In the following sections 
these two geometries will be referred to as the "lens" and 
the "pyramid" . The results of our calculations are shown 
in Table [v| and Fig. g 

A Confined electron states 

Figure ^ shows the calculated square of the envelope 
function for the electron states in the pyramidal and lens 
shaped InAs/GaAs quantum dots. For the lens shaped 
dot, the electron states can be approximately interpreted 
as eigenstates of the L z operator Q. Here we plot only 
the first 6 bound states corresponding to l z = 0, ±1 and 
±2. The first state en, has l z — and is commonly de- 
scribed as s-like as it has no nodes. The e\ and ei states 
have l z = ±1, and arep-like with nodal planes (110) and 
(110). The e3, e± and es states have l z = ±2 and respec- 
tively and are commonly described as d x 2_ y 2, d xy and 2s 
respectively. Due to the underlying zincblende atomistic 
structure, the Coo symmetry is reduced to Ci v . Hence, 
the en to es states correspond to the ax, b\, a%, ai and 
ai irreducible representations of the Ci v group, rather 
than eigenstates of L z . This allows states en,e3 and es 
to couple. This coupling is evident, for example, in the 



larger charge density along [110] compared to [110] in the 
e-i state, due to its coupling with e\. The observable ef- 
fect of this Ci v symmetry is to split the e\ and ei p-states, 
S P p, and the e% and e^ d-states, Sdd- The alignment of the 
ei and ei p-states states along the [110] and [110] direc- 
tions also results from the underlying zincblende lattice 
structure. Note, this analysis neglects the effects of the 
spin-orbit interaction which reduces the Ci v group to a 
double group with the same single representation for all 
the states. In our calculations the spin-orbit interaction 
is included, but is produces no significant effects for the 
electron states. 

The electron states in the pyramidal dot also belong 
to the Civ group and show a one-to-one correspondence 
with those in the lens shaped dot. However, there are 
only 5 bound states in the pyramidal dot due to its 
smaller size. Here we define an electron state as bound 
if its energy is below that of the unstrained, bulk GaAs 
conduction band edge. 

The calculated values of the s-p and p-d energy spac- 
ings, 5 sp , and, 5 p d, for the lens and pyramidal shaped 
dots are 65 and 68 meV and 108 and 64 meV respec- 
tively. The splitting of the two p states, 6 PP = ei — e\ 
are 2 and 26 meV respectively. The calculated values 
of the electron binding energy, AE(e), are 271 and 171 
meV respectively. The electron-electron direct Coulomb 
energies, J|g eo , Jlf^ and J| e ei in the lens and pyrami- 
dal dots are calculated as 32, 25 and 25 meV and 40, 35 
and 36 respectively. On applying a magnetic field in the 
growth direction, we calculate an increase in the splitting 
of the two p states (ei — ex) in the lens shaped dot from 2 
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to 20 meV. Details of this magnetic field calculation will 
be given in a future publication^]. Finally, the energy 

of the electron wetting layer level, AEyy L , with thick- 
nesses of 1 and 2 ML is 15 and 24 meV below the CBM 
of unstrained bulk GaAs. 



B Confined hole states 

Figure [| shows calculated wavefunctions squared for 
the hole states in pyramidal and lens shaped InAs/GaAs 
quantum dots. Unlike the electron states, the hole states 
cannot be approximated by the solutions of a single band 
Hamiltonian. Instead there is a strong mixing between 
the original bulk Bloch states with T Sv and sym- 
metry. The larger effective mass for holes results in a re- 
duced quantum confinement of the hole states and conse- 
quently many more bound hole states. Only the 6 bound 
hole states with the highest energy are shown in Figure ||. 

The calculated values of the hg-hi , hi -hi and hi-h-$ 
hole level spacings for the pyramidal and lens shaped dots 
are 8,7 and 6 meV and 15, 20 and 1 meV respectively. 
The calculated hole binding energies, AE(e), are 194 and 
198 meV. We calculate the highest energy hole level in 
pure In As wetting layers, AE^ L , with thicknesses of 1 
and 2 ML to reside 30 and 50 meV above the VBM of 
unstrained bulk GaAs. The hole-hole Coulomb energies, 
jf£ ho , are 25 and 31 meV. 



C Electron-hole excitonic recombination 

Figure]^ shows our calculated single exciton absorption 
spectrum for a pure InAs, lens shaped dot with a base of 
252 A and a height of 35 A. The energies of each of the ab- 
sorption peaks are calculated from Eq. (|l4|) . The ratios of 
the dipole matrix elements for light polarized along [110] 
and [110] are calculated from Eq . (|l6|) . Figure Q illus- 
trates that, for a lens shaped dot, both the conventional 
e; — * hi transitions and additional, ex — hi, e.% — hi, e^ — hi 
and ex — /13 transitions are strongly allowed. The ratio of 
the polarization anisotropies, A, are shown in Tabic IV. 



As a result of the circular symmetry of the lens shaped 
dot, we calculate a polarization ratio of A = 1.03 for the 
eo — ho transition. This value is in contrast to that cal- 
culated value for a pyramidal dot of A = 1.2[^3). For the 
higher angular momentum transitions we find larger de- 
viations from unity. The magnitude of the ratios, follows 
the polarization of the wavefunctions shown in Fig. ^. 
For example we find ratios greater and then less than 
unity for the ex — hi and ei — hi transitions, as reflected 
by the elongations of the ex, hi and e%, hi wavefunctions 
along the [110] and [110] directions. 

We calculate ground state electron-hole direct 
Coulomb energies, Jf^^, of 37 and 31 meV in the lens 
shaped and pyramidal dots. The calculated ground state 
electron- hole exchange energies, K^ h are an order of 
magnitude smaller, with values of 3 and 0.2 meV. These 




FIG. 4: Top view of the calculated electron and hole 
wavefunctions squared for lens and pyramidal shaped 
InAs quantum dots embedded in GaAs, with bases of 252 
and 113 A and heights of 25 and 56 A. The yellow and 
green isosurfaces represent 20 and 60 % of the maximum 
charge density. 



yield excitonic band gaps of 1.03 and 1.12 respecti vely . 
The calculated polarization anisotropy ratios [Eq.(fL6[)] 
for light polarized along [110] and [110] directions are 
A = 1.03 and 1.20 for the lens and pyramidal shapes 
respectively. The calculated excitonic dipoles [Eq.(|i7|)] 
are -3.1 and 0.16A respectively. A positive dipole is de- 
fined as the center of the hole wavefunction being located 
above the center of the electron wavefunction. 



VI. ANALYSIS OF PERTINENT 
EXPERIMENTAL MEASUREMENTS 

A The intra-band s-p and p-d electron energy 
spacings 

Measurements of the spacing between the eo and ex 
like electron levels (s-like and p-like) are based on in- 
fra red absorption. For the lens shaped dots, Fricke et. 
aZ.Q] load electrons into the dots by growing a sample 
consisting of a n-type doped layer, a tunneling barrier, 



TABLE IV: Calculated polarization anisotropy, A 
■Pi 10 : P1T0' f° r l ens sna P e d and pyramidal Ga^Ini-; 
quantum dots embedded within GaAs. 



Fig 5. Williamson ef. al. 
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Geometry 
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en — ho 


1.03 


1.2C 


ei - hi 


0.82 


2.4C 


ei — hi 


1.27 


0.52 


e 3 - h a 


0.73 


4.2f 


B4 — fl4 


1.23 


0.65 



a layer of InAs/GaAs lens shaped dots, a GaAs sp< 
and a GaAs/AlAs short period superlattice (SPS). 
applying a voltage between the n-doped layer at the I 
torn of the sample and a Cr contact grown on tor 
the SPS, electrons are attracted from the n-doped Le 
into the InAs dots. Infra-red photons were used to 
cite electrons from the occupied eo level into the ei le 
Neglecting the small exchange energy, the energy dii 
ences for the e\ — eo excitations when 1 and 2 electr 
are present in the dot are 



E 01 [e\]-E al [el\ = {e ei 
Eoi[eo<z\] - E 02 [el] = (e ei 



r jee _ jee ] 
L^ei.eo e ,e J 



The first of these energy differences yields a direct n 
surement of the s — p energy spacing, 8 sp , of 49.1 m 
The second energy difference was measured at 50.1 m 
Drexler et. al.B also used infra red transmission sj 
troscopy to measure the energy spacing, S sp = 41 m 
Pan et. al. |2l], Q have also performed infra-red abs( 
tion measurements on truncated pyramidal dots wit 
base of 180 A and height of ~60 A. In these experi- 
ments, no gate voltage is applied, and therefore the exci- 
tations take place from the ground state of the samples, 
Eqo. They observe multiple infra-red absorption peaks 
between 89 and 103 meV. These could be associated ei- 
ther with the s-p spacing of the electron levels or spacings 
of the hole states (see below) . 

Itskevich et. aZ.|l9| perform high power PL measure- 
ments of pyramidal dots with a base of 150 A and a 
height of 30 A . This high power excitation is able to 
simultaneously load multiple excitons into the dots. Due 
to state filling, these multiple excitons will occupy ground 
state (eo,/io) and higher (ei, ei, ho, hi, h^) single par- 
ticle levels. Therefore the PL measurements are able 
to simultaneously measure recombination between elec- 
trons occupying the eo,ei,e2 and e^ levels with holes 
in the ho,hi,h2 and /13 levels. In general, to describe 
the total energy differences associated with decay from 
N to N — 1 excitons occupying a dot requires a treat- 
ment that includes the exchange and correlation be- 
tween multiple occupational configurations of the N and 
N — 1 excitonic states. Such a configurational interac- 
tion (CI) approach has previously been considered for 
model parabolic dots|l| and will be discussed for realistic 
dots in a future publication |53|]. For the purposes of this 
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FIG. 5: The single exciton absorption spectrum for a 
pure InAs, lens shaped dot with a base of 252 A and a 
height of 35 A. The absorption peaks are calculated from 
Eq. (pT|) . The ratios of the dipole matrix elements for 
light polarized along [110] and [110] are calculated from 
Eq.©. 



discussion, we limit ourselves to discussing the energy 
differences associated with the lowest energy configura- 
tions on the TV exciton state, i.e. those predicted by the 
Aufbau principle. Within this approximation, the peaks 
in the high power, PL spectra can be interpreted as cor- 
responding to [see Eq.(15) where exchange is neglected] 



Peak 1: Eu[h\e\] - Eqq = (e eo - e ho ) - J' 



eo,ho 



Peak 2: E 33 {h 2 h\e 2 oe\} - E 22 [h 2 e 2 ] = (e ei - e hl ) - jg M 

1 o [ jeh jeh ■ jee ■ jhh 1 

Peak 3: E 77 [/ig/if^/i^ege?e|eJ] - E^[hlh\h\e\e\e\\ 



c e3 
2 



jeh 



£h 3 ) 

rt \ ^ [ Teh jeh 

~ Z / , i^ J e 3 M l " J ei,h 3 



hh 

h 3 ,hi 



.(19) 



Note, peak 3 is not assigned to a recombination from e2 
to h 2 as this is almost degenerate with peak 2. Itske- 
vich et. al. assume that (i) the Coulomb integrals in the 
square brackets on the right hand side of Eq.([19|) cancel, 
(ii) that Jf h = Jf h = Jf h , and (hi) that the hole 
spacings, Shohi , $h 1 h 2 > are small compared to the electron 
level spacings. With these assumptions the spacings be- 
tween peaks 1 and 2 and peaks 2 and 3 can be assigned 
to the s-p and p-d energy spacings. They find spacings 
5 sp and S pe i of 75 and 80 meV respectively. Our calcula- 
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tions suggest that assumptions (i), (ii) and (hi) probably 
introduce errors of ~ ±10, ~ ±5 and +10 meV respec- 
tively. The neglect of exchange interactions in the above 
discussion also introduces an error of ~ ±5 meV. 



B The intra-band electron p level splitting 

For the lens shaped dots, the Capacitance Voltage 
spectroscopy of Fricke et. al. [j can be used to estimate 
the splitting of the p states, S PPl by loading two electrons 
into the dot and exciting them using FIR spectroscopy. 
The relevant energy differences are 

E 02 [ele\] - E 2[e 2 ] = (e ei - e eo ) + [j™ eo - J| %J 

E 02 [e l e\] ~ E 02 [e 2 ] = (e e2 - e eo ) + [J e e 2 % - J e e % ] (20) 

By assuming J* e . = Jt e „ , the difference in the two 

J o ei,eo fi2i e o' 

above expressions yields the energy spacing e e2 — e ei . 
They hnd a value of ~2 meV. To measure the effect of 
a magnetic field on the splitting of the p states Fricke 
et. use infra-red absorption to measure the above 

energy differences in an applied magnetic field. At a field 
of 15 Tesla they measure an energy spacing of 19 meV. 
More theoretically, Dekel et. a?[Q have demonstrated 
that one must assume a splitting of the p states to ex- 
plain the number of multi-exciton levels observed in their 
single dot, micro PL measurements. 

C The intra-band hole energy spacings 

There are currently no measurements available for 
the energy spacings between the hole states in the lens 
shaped dots. Itskevich et. aL[jl9j have performed high 
power PL measurements of dots estimated to be square 
based, truncated pyramids with a base of 150A and a 
height of 30A , under hydrostatic pressure to estimate 
hole level spacings. At hydrostatic pressures above 55 
kbar, they measure the PL associated with transitions 
from the Ai c -state in the GaAs matrix to the ho and h\ 
levels in the dots. They estimate a hole level spacing of 
~15 meV. However, the nature of the initial and final 
hole states is unclear. 

Sauvage et. aL[|22| use polarized photoinduced in- 
traband absorption spectroscopy to measure the energy 
spacing between the lowest hole state, hooo: and the hole 
state with a single node in the growth direction, hooi- 
This corresponds to [see Eq.(|l5|)] 

E\i[hl m el] - i?n[ft-ooi e o] = (e^ooi ~ e 'iooo) 

By assuming that Jf 1 h — jf 1 h Sauvage et. al. esti- 
mate the /iooi~^ooo spacing to be ^120 meV. Note h oi is 
almost certainly higher in energy than states with nodes 
perpendicular to the growth direction (/iqio and /iioo) due 



to the smaller dimension of the dot in the growth direc- 
tion. Consequently, the energy difference hooi — ^ooo is 
not the spacingof the first two hole states ho-hi. 

Tang et. al. [p4| measure activation energies for excita- 
tions from ho and hi to the hole wetting layer of 48 and 
30 meV respectively, implying an ho~h\ spacing of ~18 
meV. 



D The electron and hole binding energies, AE(e) 
and AE(h) 

There have been no direct measurements of the elec- 
tron or hole binding energy for lens shaped InAs dots. 
However, it has been measured in other dots by sev- 
eral groups using a range of techniques. Berryman et. 
al.U% placed pyramidal InAs dots estimated to have a 
base of 100 A and height of 15 A in a p-n junction and 
measured the temperature dependence of the ac conduc- 
tance as a function of frequency. These measurements 
predict a hole binding, AE(h), energy of ^240 meV. 
When subtracted from the bulk GaAs band gap, this 
yields a value for the electron binding energy, AE(e), of 
~60 meV. The authors obtain similar results from tem- 
perature dependent Hall measurements of thermal hole 
trapping. Itskevich et. ffl^ flSt measured the pressure at 
which PL measurements could detect a T-X crossing in 
pyramidal InAs / GaAs quantum dot samples. By extract- 
ing these PL measurements back to zero pressure they 
were able to extrapolate a value for the electron bind- 
ing energy, AE(e), of ^50 meV. Itskevich et. aZ.jl9) also 
used high pressure PL to measure the energy difference 
between the X\ c level in bulk GaAs and the ho level in 
the quantum dots. By extrapolating this value back to 
zero pressure, they predict a value for the hole binding 
energy, AE(h), of ~250 meV. Brunkov et. aZ. [^3j per- 
formed capacitance-voltage spectroscopy measurements, 
which when fitted to a capacitance model for their dot ge- 
ometry predict an electron binding energy of 80 meV for 
dots with bases of 250 A. Tang et. aL|5J] measured the 
spacing of the electron wetting layer to both the GaAs 
CBM and the eo level and hence deduced a value for the 
electron binding energy, AE(e), of ~80 meV, for dots 
with an estimated base of 130 to 170 A. 



E The position of the electron and hole wetting 
layer level 

The presence of a distinct wetting layer signal in the 
PL spectra of a sample of self assembled quantum dots 
is the hallmark of a high quality sample. In samples 
where the wetting layer has "dissolved" due to the growth 
conditions, it is likely that the geometry and composition 
of the quantum dots will also have dramatically altered 
from their uncapped state. 

In the lens shaped InAs dots Schmidt et. qL[|Tcj] ob- 
serve PL emission from the ground state of the wetting 



11 



layer at 1.34 eV. Photovoltage measurements M] on the 
same samples show a strong peak corresponding to ab- 
sorption into the ground state of the wetting layer at 1.35 
eV. There are currently no measurements available for 
the position of the individual electron and hole wetting 
layers in the lens shaped InAs/GaAs quantum dot sam- 
ples. In Ref. |^3| Sauvage et. al. grow lens shaped InAs 
dots with an estimated base of 150 A and a height of 30 
A on a substrate that is n-doped with silicon. This n- 
doping loads electrons into the eo state in the dot, which 
they excite into the wetting layer using infra-red excita- 
tion. In these samples they estimate the wetting layer 
to be 150 meV above the eo level. In Ref. 22 Sauvage 
et. al. load electrons into the eo state of similar dots 
using an optical interband pump. Subsequent infra-red 
absorption places the wetting layer 190 meV above the 
eo state. Tang et. al. 54 measure thermal transfer of 
holes to the wetting layer, and obtain a spacing between 
the ho level and the hole wetting layer, AE^) L of ^48 
meV. They also measure thermal transfer from an ex- 
cited state, possibly involving hi, which places the hole 
wetting layer ~30 meV below the hi level. 



F The number of confined electron and hole states 

The actual number of confined electron and hole states 
in a self assembled InAs/GaAs quantum dot depends on 
the size and composition of the dot. Early single band, 
effective mass calculations^] for pure InAs pyramidal 
dots with a base of 120 A and height 60 A predicted only a 
single bound electron state and several bound hole states. 
Consequently, many experiments were then interpreted 
in this light. More accurate multi-band /c.p]55|, |5(| O] 
and pseudopotential|35[ |36| calculations have predicted 
5 or more bound electron states in the same dots. 

The high power PL experiments of Itskevich et. al. |l^] 
show the gradual disappearance of 5 peaks as a func- 
tion of external pressure. This is interpreted as direct 
evidence for 5 confined electron levels in their samples. 
The single dot, multi-exciton measurements of Dekel et. 
a/, [hi require the assumption of at least 3 bound elec- 
tron states to explain their experimental spectra. Sim- 
ilarly, the capacitance-voltage spectroscopy of Fricke et. 
aZ.j?]] shows two peaks corresponding to the capacitance 
of s-like states and the nearly degenerate p-like states, 
providing evidence for at least 3 bound electron states. 



Electron and hole Coulomb and Exchange 
Interactions 



By loading multiple electrons and holes into quantum 
dots it is possible to measure the Coulomb and exchange 
interactions between these additional electrons and holes. 
The magnitude of these interactions is a function of the 
shape of the electronic wavefunctions (see section ||) and 



provides an additional quantity to test the accuracy of 
theoretical models. 

To study electron-electron interactions, Fricke et. 
al.f^ use the same experimental setup discussed in sec- 



tion VIA. From Eq.(|l5[) we see the energy differences 
corresponding to the peaks in the Capacitance Voltage 
(CV) spectra associated with loading one and two elec- 
trons into the eo level in the dots is 



E02 H] 



En 



Eqi 



Ji, 



(22) 



The electron-electron Coulomb interaction, J^ eg , can 
therefore be directly measured as the splitting of these 
two CV peaks. They find a value of e ~23 meV. 
From Eq.(|l8|) we see that 



Tee Tee 

J e ,ei ~ J e ,e 



(50.1 - 49.1) = 24meV 



(23) 



Finally by fitting 4 equidistant bell curves to the CV 
spectra corresponding to loading 3,4,5 and 6 electrons 
into the dots, an approximate value for the charging 
energy between the p states, , of ~18 meV is ob- 

tained. From Eq.(|i"5j) we that the spacing Eqa — E 03 = 
JZU + 2J£ ei1 while E 06 - E 05 = 3J e e f iei + 2J- ej and 
hence the approximation of equidistant peaks will intro- 
duce some error into this estimate for Jf e „ . 



H Electron-hole excitonic recombination 

To our knowledge there have so far been no measure- 
ments of the polarization anisotropy in the lens shaped 
dots discussed here. The polarization anisotropy for 
e, — hi excitonic recombination in InAs/GaAs was mea- 
sured by Yang et. a/.|l^, [l(| , who find a ratio of 
^e ,h = 1-2 and X e2 ,h 2 = 1-3 for InAs dots whose ge- 
ometry is measured to be formed by four {136} faceted 
planes with bases ranging from 150 to 250 A and a base 
to height ratio of 4:1. Yang et. al.^^ [ll| have performed 
k.p calculations for this dot geometry, which include the 
"geometric factor" but not the "atomic symmetry fac- 
tor" discussed in section IV. They find A e „,h = 1.8 and 
\ e2 ,h 2 = 3.5. The authors suggest the k.p simulations of 
the measured polarization ratio can be used to deduce 
the geometric shape anisotropy. However, we demon- 
strate here that when the "atomic symmetry factor" is 
included an anisotropy of X eo ,h = 1.3 is obtained even 
for a square based pyramid. Thus k.p simulations lacking 
the "atomic symmetry" factor can not be used to reliably 
deduce the geometric shape anisotropy. 

To our knowledge there have so far been no measure- 
ments of the excitonic dipole in lens shaped InAs dots. 
Fry et. al. [ pOfl used photocurrent spectroscopy within an 
applied electric field to measure the excitonic dipole mo- 
ment, dh it ej, [Eq.(|l7|)]. They find the center of the hole 
wavefunctions to be located ~4 A above the center of the 
electron wavefunction (positive dipole). Fry et. al. 63] 
also perform single-band, effective mass calculations, in 
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an attempt to isolate the origin of this dipole. They 
predict that in the absence of alloying the dipole is - 
3 A, i.e. the opposite sign , but a linear composition 
profile with Gao.5Ino.5As at the base and pure InAs at 
the top of a truncated pyramid with a base of 155 A, 
and height 55 A, reproduces the correct dipole of 4 A. 
They suggest that this alloying profile explains the ob- 
served dipole. We have repeated these calculations and 
confirm that, within a single-band model, such a com- 
position profile causes both electrons and holes to move 
up in the dot compared to their positions in a pure InAs 
dot. The heavier effective mass of the holes, results in 
less kinetic energy associated with confinement at the 
top of the dot and hence the holes move up more than 
electrons on the introduction of Ga, producing the cor- 
rect dipole. However, when we repeat these calculations 
in the more sophisticated, multi-band LCBB basis used 
here, we find significant heavy hole-light hole mixing in 
the ho state, which acts to reduce the above effect and 
produce a smaller dipole of — 1 A, in contradiction with 
experiment (+4 A). We therefore conclude that in a more 
complete calculation the shape and alloy profile postu- 
lated by Fry et. al.^H^ does not produce the observed 
excitonic dipole. 



VII. COMPARISON OF EXPERIMENT AND 
THEORY 

In Table [v| we show the results of our calculations for 
pure InAs, lens shaped quantum dots embedded within 
GaAs [column (a)]. Table [v] also shows the experi- 
mentally measured splittings of the electron levels, the 
electron-electron and electron-hole Coulomb energies, the 
magnetic field dependence and the excitonic band gap 
measured in Refs.jTj and |ll|. The agreement between 
the measured energy level spacings, Coulomb energies 
and magnetic field response with our theoretical lens 
shaped model is generally good. Both the model and 
experiment find (i) a large spacing, 8 sp , (—50-60 meV) 
between the s-like eo state and the p-like e\ state, (ii) a 
small spacing, <5 PP , (—3 meV) between the two p-like e\ 
and ei states and (iii) a large spacing (—55 meV) between 
the p-like e-i state and the d-like state. 

These electron level spacings are similar to those found 
for pyramidal quantum dots[p5| (see Table [y]). However, 
in the pyramidal dot, the spacings of the two p-like and 
d-like states, d pp , ddd, is larger (26 and 23 meV) as a result 
of the lower C2„ symmetry of a zincblende pyramid. Both 
the model and experiment also find similar values for the 
Coulomb energies, J(eoeo) and J(eoho) (—25 meV). 

The calculated hole binding energy of AE(h) — 193 
meV is in good agreement with those of Berryman et. 
al.$\J§ (-240 meV) and Itskevich et. a/.Q (-250 meV). 
Our calculated electron binding energies, AE{h), are con- 
siderable larger (271 meV) than those of Berryman et. 
ai.@ (-60 meV) and Tang et. a/.|| (-80 meV). We 
attribute this difference to the larger size of our dots. 



The assumption of a pure InAs dot also affects the com- 
parison. The agreement improves when we include Ga in 
our dots (see section VII B). 

The calculated electron-electron and electron-hole 
Coulomb energies are in reasonable agreement with those 
extracted from Refs.Q and [[fl|. For the integrals 

J !oe > J e e ei > J eTei and J toh we calculate values of 31, 25, 
25 and 37 respectively, compared to measured values of 
23, 24, 18 and 33.3 meV 

The calculated polarization anisotropics, A, for the 
eo — ho recombination in lens and pyramidal shaped, pure 
InAs dots are A = 1.03 and 1.2 respectively. A future 
measurement of this anisotropy ratio for lens shaped dots 
would provide an important piece of evidence for deter- 
mining the detailed geometry of the dots. 

In the lens shaped dot we find a difference in the av- 
erage positions of the ho and eo states, dhi, ej , of around 
1 A. This is smaller than the value we calculate for a 
pyramidal quantum dot, where we find the hole approx- 
imately 3.1 A higher than the electron. 

In summary, the assumed lens shaped geometry, with a 
pure InAs composition produces a good agreement with 
measured level splitting, Coulomb energies and magnetic 
field dependence. A closer inspection of the remaining 
differences reveals that the calculations systematically 
overestimate the splittings between the single particle 
electron levels (S sp : 65 vs. 50 meV, <5 p( j:68 vs. 48 meV) 
and underestimates the excitonic band gap (1032 vs. 1098 
meV). 



A Pure InAs dots: The effects of shape and size 

Focusing on the lens shape only, we examine the effect 
of changing the height and base of the assumed geome- 
try. Calculations were performed on similar lens shaped, 
pure InAs dots where (i) the base of the dot was increased 
from 252 to 275A, while keeping the height fixed at 35A, 
[column (b)] and (ii) the height of the dot was decreased 
from 35 to 25A, while keeping the base fixed at 252A, [col- 
umn (c)]. It shows that decreasing the height of the dot 
increases the quantum confinement and hence increases 
the splittings of the electron and hole levels (S sp : from 65 
to 69 meV and 8h ,h-i ■ from 8 to 16 meV). Decreasing the 
height of the dot also acts to increase the excitonic band 
gap from 1032 to 1131 meV by pushing up the energy 
of the electron levels and pushing down the hole levels. 
Conversely, increasing the base of the dot decreases both 
the splittings of the single particle levels (S sp : from 66 
to 61 meV) and the band gap (1032 to 1016 eV). These 
small changes in the geometry of the lens shaped dot have 
only a small effect on electronic properties that depend 
on the shape of the wavefunctions. The electron-electron 
and electron-hole Coulomb energies remain relatively un- 
changed, the magnetic field induced splitting remain at 
20 meV, the polarization anisotropy, A, remains close to 
1.0 and the excitonic dipole, dh it e remains negligible. 
In summary, reducing either the height or the base of 
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the dot increases quantum confinement effects and hence 
increases energy spacings and band gaps, while not sig- 
nificantly effecting the shape wavefunctions. 



B Interdiffused In(Ga) As/GaAs lens shaped dots 

We next investigate the effect of changing the compo- 
sition of the quantum dots, while keeping the geometry 
fixed. There have recently been several experiments [20, 
|25| , |26| suggesting that a significant amount of Ga dif- 
fuses into the nominally pure InAs quantum dots during 
the growth process. We investigate two possible mecha- 
nisms for this Ga in-diffusion; (i) Ga diffuses into the dots 
during the growth process from all directions producing 
a dot with a uniform Ga composition Ga^Ini-^As, and 
(ii) Ga diffuses up from the substrate, as suggested in 
Ref. j2(|. To investigate the effects of these two meth- 
ods of Ga in-diffusion on the electronic structure of the 
dots, we compare pure InAs dots embedded in GaAs 
with Ga^Ini-^As, random alloy dots embedded in GaAs, 
where the Ga composition, x, (i) is fixed at 0.15, [column 
(d)] and (ii) varies linearly from 0.3 at the base to at 
the top of the dot, [column (e)]. 

The electronic structure of these dots is compared in 
Table [v|. ft shows that increasing the amount of Ga in 
the dots acts to decrease the electron level spacings (S sp : 
from 65 to 58 and 64 meV for x — 0.15 and x = 0.3 to 
x = respectively). It also acts to increase the excitonic 
band gap from 1032 to 1080 and 1125 meV respectively. 
The electron binding energy, AE(e), is decreased by the 
in diffusion of Ga (from 271 to 209 and 192 meV), while 
the hole binding energy, AE(h), is relatively unaffected. 
This significant decrease in the electron binding energy 
considerably improves the agreement with experiments 
on other dot geometries [p7| [5^ ]. 

As with changing the size of the dots, we find that Ga 
in-diffusion has only a small effects on properties that de- 
pend on the shape of the wavefunctions. The calculated 
electron-electron and electron-hole Coulomb energies are 
almost unchanged, while the average separation of the 
electron and hole, dh i: ej, increases from 0.16 to to 0.5 
and 1.2 A and the polarization ratio, A, and magnetic 
field response are also unchanged. 

Table shows that the dominant contribution to the 
increase in the excitonic band gap and reduction in elec- 
tron binding energy, results mostly from an increase in 
the energy of the electron levels as the Ga composition 
is increased. This can be understood by considering the 
electronic properties of the bulk Ga 2 .In 1 _ a; As random al- 
loy. The unstrained valence band offset between GaAs 
and InAs is ~ 50 meVp7|, while the conduction band 
offset in ~ 1100 meV and hence changing the Ga compo- 
sition, x, has a large effect on the energy of the electron 
states and only a small effect on the hole states. In sum- 
mary, the effect of Ga in-diffusion is to reduce the spacing 
of the electron levels while significantly increasing their 
energy and hence increasing the band gap. We find that 



only the average Ga composition in the dots is important 
to their electronic properties. Whether this Ga is uni- 
formly or linearly distributed throughout the dots has a 
negligible effect. Note, in Ref.[^o) it is suggested that a 
linear composition profile is required to produce an ex- 
citonic dipolc moment in agreement with that measured 
by stark experiments. For the lens shaped geometry dis- 
cussed here there have so far been no such measurements 
of the dipole, but our calculations suggest that it should 
be small (~lA). 



VIII. DISCUSSION 

The effects of changing the geometry of the lens shaped, 
pure InAs dots on the single particle energy levels can be 
qualitatively understood from single band, effective mass 
arguments. These predict that decreasing any dimension 
of the dot, increases the quantum confinement and hence 
the energy level spacings and the single particle band 
gap will increase. Note that as the dominant quantum 
confinement in these systems arises from the vertical con- 
finement of the electron and hole wavefunctions, changing 
the height has a stronger effect of the energy levels than 
changing the base. In this case decreasing the height by 
10 A has a much stronger effect on the energy spacings 
and on the band gap than increasing the base by 23A. 

As increasing(decreasing) the dimensions of the dot 
acts to decrease(increase) both the level spacings and the 
gap, it is clear that changing the dot geometry alone will 
not significantly improve the agreement with experiment 
as this requires a simultaneous decrease in the energy 
level splittings and increase in the band gap. However, 
Ga in-diffusion into the dots acts to increase the band gap 
of the dot while decreasing the energy level spacings. Ta- 
ble |V| shows that adopting a geometry with a base of 275 
A and a height of 35 Aand a uniform Ga composition of 
Gao.15Ino.85As produces the best fit to the measurements 
in Refs.§ and @. 

In conclusion, our results strongly suggest that to ob- 
tain very accurate agreement between theoretical models 
and experimental measurements for lens shaped quan- 
tum dots, one needs to adopt a model of the quantum 
dot that includes some Ga in-diffusion within the quan- 
tum dot. When Ga in-diffusion is included, we obtain an 
excellent agreement between state of the art multi-band 
pseudopotential calculations and experiments for a wide 
range of electronic properties. We are able to predict 
most observable properties to an accuracy of ±5 meV, 
which is sufficient to make predictions of both the geom- 
etry and composition of the dot samples. 
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TABLE V: Calculated single particle electron and hole energy level spacings, electron and hole binding energies, 
AE(e, h), electron-electron and electron-hole Coulomb energies, excitonic band gap all in meV, exciton dipole moment 
and polarization anisotropy for lens shaped and pyramidal Ga^Ini-^As quantum dots embedded within GaAs. 
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